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ABSTRACT 

We have carried out JHK polarimetric observations of eleven dusty young stars, 
by using the polarimeter module IRPOL2 with the near-IR camera UIST on the 3.8-m 
United Kingdom Infrared Telescope (UKIRT). Our sample targeted systems for which 
UKIRT-resolvable discs had been predicted by model fits to their spectral energy 
distributions. Our observations have confirmed the presence of extended polarized 
emission around TW Hya and around HD 169142. HD 150193 and HD 142666 show the 
largest polarization values among our sample, but no extended structure was resolved. 
By combining our observations with HST coronographic data from the literature, we 
derive the J- and H-band intrinsic polarization radial dependences of TW Hya's disc. 
We find the disc's polarizing efficiency is higher at H than at J, and we confirm that 
the J- and H-band percentage polarizations are reasonably constant with radius in 
the region between 0.9 and 1.3 arcseconds from the star. We find that the objects for 
which we have detected extended polarizations are those for which previous modelling 
has suggested the presence of flared discs, which are predicted to be brighter than flat 
discs and thus would be easier to detect polarimetrically. 

Key words: stars : evolution - Circumstellar matter - planetary systems: protoplan- 
etary discs - infrared : stars - techniques: polarimetric 



1 INTRODUCTION 

The discovery of the Vega- phenomenon in 1984 is consid- 
ered one of the most important achievements of the IRAS 
space telescope. It discovered that several main sequence 
(MS) stars, including Vega, exhibited large mid- and far- 
infrared excesses that could not be ascribed to pure photo- 
spheric emission (Aumann et al. 1984). Such excesses were 
attributed to a disc or ring of solid particles surrounding 
the stars, later termed 'debris-discs', and great interest was 
taken in these objects because of their relevance to the for- 
mation of planetary systems. 

The case of young low- and intermediate-mass sys- 
tems, such as T Tauri and Herbig Ae/Be (HAeBe) stars, 
is of particular interest (Yudin, 2000; Waters & Waelkens, 
1998). These stars are considered to be precursors of classic 
Vega-like stars, such as (3 Pictoris, e Eridani, HR 4796A, Fo- 
malhaut and Vega itself. While planetary bodies are thought 
to have already formed in Vega-like systems (Aumann et al. 
1984; Beust et al., 1996; Natta et al., 2004), the circum- 
stellar environments of T Tauri and HAeBe stars probably 
represent an early phase of planet formation (Bouwman et 
al., 2001). 

From high resolution optical spectra of a sample of 
young dusty stars, the majority in common with our cur- 



rent sample and classified as HAeBe or T Tauri stars on the 
basis of their emission line spectra, Dunkin, Barlow & Ryan 
(1997a; 1997b) found that analysis of their absorption line 
spectra yielded surface gravities typical of MS stars, suggest- 
ing that they had already arrived on the MS. These type of 
objects have been referred to as 'transition' objects (Mal- 
fait, Bogaert & Waelkens, 1998) or 'dusty' Vega-excess stars 
(Dent, Greaves & Coulson, 2005) and they are nowadays 
thought to play a pivotal role in our understanding on the 
transition between the pre-MS (T Tauri and HAeBe stars) 
and the MS (Vega- like stars). 

Many issues remain unclear regarding the circumstel- 
lar structures of these 'transition' objects, in particular, 
whether they are surrounded by flat discs or spherical en- 
velopes is still not determined (Vink et al., 2002). Based 
on their spectral energy distributions, Meeus et al. (2001) 
introduced a classification of HAeBe stars into two groups: 
Group I objects have an almost flat or even rising 20-100 /im 
spectral energy distribution (SED), which can be well mod- 
elled by a flaring disc (Dullemond & Dominik, 2004) . On the 
other hand, Group II objects have a much bluer IR SED, and 
are believed to be associated with self-shadowed non-flared 
discs, which are predicted to be smaller in spatial extent 
than those of Group I sources. 

High-contrast, high-spatial resolution, imaging close to 
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bright central stars is required to test the validity of these 
models. Despite the efforts of many groups, limitations on 
instrumental capabilities have to date kept the number of 
directly imaged discs to a very low- number (see Zuckerman 
(2001) for a review). When only SEDs are available, not all 
parameters can be uniquely determined, e.g. the grain-size 
distribution and the spatial density profile are usually not 
separable. 

Imaging polarimetry provides a powerful technique for 
detecting dust-discs around bright stars. Since only the light 
from the discs is expected to be polarized, the bright central 
stars are automatically suppressed in polarized light images. 
The utility of imaging polarimetry was demonstrated over 
ten years ago when Gledhill, Scarrott & Wolstencroft (1991) 
detected polarized emission from the nearby Beta Pictoris 
disc (> 1000 AU in radius at a distance of 19.3 pc; Crifo 
ct al (1997)). The light was found to be 17% polarized and 
extended up to 30 arcseconds from the star (R- waveband). 
Kuhn, Potter & Parise (2001) also used imaging polarimetry 
to detect polarized light around the A5Ve star HD 169142. 
More recently, using the Very Large Telescope (VLT), Apai 
et al. (2004) tested the power of imaging polarimetry on an 
8-m class telescope, detecting the presence of polarized emis- 
sion as close as 0.1 arcsecond (~ 6 AU) from the classical T 
Tauri star TW Hya. 

We have carried out a JHK imaging polarimetry sur- 
vey for a sample of ten 'transition' objects (so called 'dusty' 
Vega-excess stars), selected from the surveys of Sylvester 
et al. (1996) and Mannings & Barlow (1998), along with 
TW Hya. Modelling of their SEDs had predicted angular 
sizes that could be resolvable at near-infrared (near-IR) 
wavelengths using the near-IR imaging instrument UIST on 
the 3. 8-m United Kingdom Infrared Telescope (Sylvester & 
Skinner, 1996; Sylvester, Skinner & Barlow, 1997). Our pri- 
mary aim is to detect the signature of extended polarized 
emission from a scattering disc around the target stars. Note 
that even if the disc is too small to be resolved with UKIRT, 
a net polarization can be observed, indicating the presence 
and orientation of a disc on the sky. 

Section 1.1 describes our target sample and basic selec- 
tion criteria. §2 reports on the observations and §3 on the 
data reduction; Results are presented in §4, and discussed 
in §5. 

1.1 Target sample 

We have selected a sample of ten young stars with dust discs 
from the surveys of Sylvester et al. (1996) and Mannings & 
Barlow (1998). All targets show mid- and far-IR excesses due 
to warm and cold dust grains. Most of them have been thor- 
oughly studied using spectrophotometric techniques. Spec- 
tral types, effective temperatures, surface gravities, elemen- 
tal abundances and overall line-emission characteristics are 
available for the majority of the sample (Dunkin, Barlow & 
Ryan, 1997a; Dunkin, Barlow & Ryan, 1997b). A summary 
of our sample is shown in Table 1, where we have listed 
respectively their names, spectral classification, infrared ex- 
cess fractions Lir/L*, distances and V-band magnitudes (as 
listed in the SIMBAD database). 

The IR excesses of our targets are all much larger 
than those of true Vega-like stars ( which are typically 
Lir/L* < 10~ 3 for Vega-like stars; Sylvester et al. (1996)). 



Their surface gravities indicate they have recently reached 
the MS (Dunkin, Barlow & Ryan, 1997a; Dunkin, Barlow 
& Ryan, 1997b). HD 123160, HD 123356, HD 141569 and 
HD 145263 have relatively low IR excesses compared to 
the rest of our sample. The remaining stars have consid- 
erably larger IR excesses (Lir/L* > 0.1), suggesting that 
they might be in a somewhat younger evolutionary state. 
HD 142666, HD 143006, HD 144432 and HD 169142 have 
near-IR excesses attributable to hot dust emission, result- 
ing in (J-H) and (H-K) colours coincident with the locus 
of HAeBe and T Tauri stars (Sylvester et al., 1996). It has 
been proposed that the dust responsible for the near-IR ex- 
cesses is a main agent for the optical polarization, and thus 
that Vega-like 'prototypes' showing almost no near-IR ex- 
cess should present low polarization values (Yudin, 2000). 
HD 142666, HD 143006, HD 144432 and HD 169142 have 
been detected in CO in the survey by Dent et al. (2005), 
indicating they still have substantial gas content. 



2 OBSERVATIONS 

The J-band data presented in this work was acquired on 
April 29th 2003, using the dual-beam polarimeter module 
IRPOL2 in conjunction with the near-IR camera UIST on 
UKIRT at Mauna Kea observatory. The resulting field of 
view was of 15x60 arcseconds, with 0.12 arcseconds per 
pixel. H- and K-band data for TW Hya, HD 150193 and 
HD 141569 were also acquired, on April 30th 2003, using 
the same observing procedure. 

Dual-beam imaging polarimetry was obtained using a 
lithium niobate Wollaston prism, along with a rotating half- 
wave plate. By observing at four different half-wave plate 
positions (separated by 22.5 degrees), linear polarimetry is 
acquired. A detailed description of this standard dual-beam 
polarimetry configuration can be found in Berry & Gledhill 
(2003). Basically, two orthogonally polarized images are si- 
multaneously recorded. The subtraction of these two images 
leads to the suppression of unpolarized light (mainly from 
the central star), and to the highlighting of polarized light 
(e.g. scattered light from the disc). 

We observed a total of 11 objects. Since most of our tar- 
gets are bright (K< 8), several short exposures (< Is) were 
co-added, in order to overcome saturation problems. Typical 
integration times ranged from 315 to 475 seconds. A non- 
destructive readout mode was used (ND+STARE). Point 
Spread Function (PSF) and flux calibrators were regularly 
interspersed, and polarized standards from the UKIRT cal- 
ibration catalogue were also observed. In order to minimise 
the effects of pixel to pixel variations in the detector on de- 
rived polarization values, each measurement was repeated 
at 3 different array (jitter) positions. 



3 DATA REDUCTION 

All data were identically reduced using the STARLINK 
data-reduction pipeline ORAC-DR*. Dark subtraction and 
flat-fielding are carried out, and the dual-beam imaging po- 
larimetry package POLPACK (Berry & Gledhill, 2003) was 

* http://www.oracdr.org/ 
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then used to perform component image alignment and to 
combine the data to form resultant I, Q and U Stokes im- 
ages. Propagation of the variance estimates from the raw 
data through the calculation provides errors on the final I, 
Q and U images, and on the derived polarized quantities P 
and PI (degree of polarization and polarized flux, respec- 
tively). 

The values of P presented in this work correspond to the 
degree of polarization averaged over 0.36 arcsecond square 
bins (3x3 pixels), after applying a 2 a cut. In this way, only 
measurements with an error in the degree of polarization of 
less than one-half of the polarization are considered. 

4 RESULTS 

4.1 J- band imaging polarimetry and photometry 

J- band photometry is presented in Table 2, together with 
2MASS photometry. The photometric values derived from 
our measurements are generally in good agreement with the 
2MASS values. HD 143006's J- band magnitude shows a 0.23 
magnitude deviation with respect to previous J- band pho- 
tometry, suggesting that the object has brightened since the 
2MASS measurements. HD 169142's J- band magnitude is 
also fainter than the 2MASS photometry, but is consistent 
with J- band measurements from Sylvester et al. (1996). 

Table 3 presents measured integrated polarizations for 
our targeted stars. In the case of resolved centrosymmetric 
patterns (as in the cases of TW Hya and HD 169142), the 
integrated polarization will appear less than the actual re- 
solved polarizations, since integrating the polarization vec- 
tors over all orientations can average to zero. As an example, 
TW Hya's maximum resolved polarization is 15 ± 0.5 % 
at 1.6 arcseconds from the star, while its integrated polar- 
ization is only 0.233 ± 0.025 %. The J- band instrumental 
polarization was estimated from observations of unpolarized 
standards, which showed percentage polarization levels of 
less than 0.5 percent. These values represent the real pre- 
cision of our polarimetric measurements (mostly residuals 
from the alignment procedure), and should again not be 
confused with the statistical errors on the integrated polar- 
izations quoted in Table 3. 

Targets for which we detect evidence of polarizing cir- 
cumstellar material are discussed individually below (TW 
Hya, HD 169142, HD 150193 and HD 142666). The remain- 
ing stars from our sample do not show any evidence of dusty 
polarizing discs. This suggests that their polzarization levels 
are too faint to be detected by our observations, as exem- 
plified by the non-detection of the disc around HD 141569. 
The intensity distribution of HD 141569's ring-like disc is 
known to peak at 3.3 arcseconds from the star (Augereau 
et al., 1999), a scale resolvable with UKIRT. However our 
J-band polarimetric observations were not able to detect it. 
Combining the HST J-band coronographic observations of 
HD 141569's disc (Augereau et al., 1999) and our polarimet- 
ric data, we derive upper-limits on the polarization level of 
HD 141569's disc. From this we conclude the disc is so faint 
that even if it was 100% polarized it would still fall below 
our detection limit. Augereau et al.'s (1999) HST observa- 
tions showed that the surface brightness of HD 141569's disc 



peaks at 0.199 ± 1.2 mjy arcsec -2 . This is one order of mag- 
nitude fainter than that of TW Hya's disc (5.7 ± 1.4 mjy 
arcsec -2 ; Weinberger et al. (2002)), and could explain why 
the disc around TW Hya could be seen in polarized light, 
while the one around HD 141569 was not. 

J-band imaging polarimetry of four of our eleven targets 
is presented in Figures 1, 2, 4 and 5 using the same format 
for each target. Polarization vector maps are superimposed 
upon total intensity contours. Total intensity grey-scales are 
also plotted, with the purpose of highlighting the position 
of the central source. The polarization vectors are oriented 
parallel to the E vector, with their length proportional to 
the degree of linear polarization, as indicated by the scale 
vector on each diagram. 



TW Hya: 

Our observations show a centro-symmetric polarization pat- 
tern (Figure 1, left panel), revealing the signature of a scat- 
tering disc seen near to face-on. The polarization values in- 
crease with increasing radius, up to ~ 15 % at 1.6 arcsec- 
onds, although these have larger uncertainties than those 
close in. We note the ~ 15 % percentage polarization val- 
ues we measure at 1.6 arcseconds include direct unpolarized 
flux from the star in the estimate of the disc's total scattered 
light intensity, and so represent a lower limit on the disc's 
intrinsic polarization. Please refer to Section 5 and Figure 
9, left panel, for an estimate on the intrinsic polarization of 
TW Hya's disc. We estimate the extension of the polarizing 
disc by comparing the polarized intensity (PI) radial profile 
of TW Hya to the total intensity radial profile of the PSF 
reference star HIP54690 (Figure 1, right panel). Extended 
polarized structure can be traced from ~ 0.4 arcseconds up 
to at least ~ 1.5 arcseconds. This is comparable to the ex- 
tension of TW Hya's polarizing disc as seen in the Ks- band 
with the VLT (1.4 arcseconds in radius, Apai et al. (2004)). 

HD 169142: 

Our UKIRT observations just marginally resolve a centro- 
symmetric pattern around HD 169142 (Figure 2, left panel). 
The PI radial profile of HD 169142 (Figure 2, right panel) 
shows the detection of an extended polarizing disc from 
~ 0.4 to about 1.2 arcseconds radius. Inspection of the Q- 
polarization image of HD 169142 (Figure 3, left panel) re- 
veals a clear modulation with angle around the star (Figure 
3, right panel). The modulation reveals the presence of an 
extended distribution of dust, as has already been shown by 
Kuhn et al. (2001). 

HD 150193: 

Recent adaptive optics H-band coronographic observations 
of the binary system HD 150193 A-B (1.1 arcseconds separa- 
tion) at the 8.2-m Subaru telescope (Fukagawa et al., 2003) 
have revealed the presence of a scattering disc around the 
primary star, extending from the edge of the coronographic 
mask (0.5 arcseconds radius, or 50 AU) to about 1.3 arcsec- 
onds. An asymmetry of the disc in the direction of the binary 
companion was detected, suggesting that the companion has 
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a distorting effect on the disc structure. HD 150193 A is 
known to be a A2 IV Herbig Ae star (Mora et al., 2001), 
whilst HD 150193 B has recently been classified as a classi- 
cal T Tauri star (Bouvier & Corporon, 2001). 

Our J-band observations show polarization vectors aligned 
at a position angle (PA) of approximately 57 degrees, mea- 
sured from North to East (Figure 4, left panel). This sort 
of pattern is usually associated with scattering in an un- 
resolved geometry. In fact, the J-band PI radial profile of 
HD 150193 A shows no extension when compared to the to- 
tal intensity radial profile of the PSF reference star (Figure 
4, right panel). 

HD 150193 presents the highest integrated percentage po- 
larization of the sample (3.15% at PA=57 degrees), of which 
a significant part is thought to be interstellar in nature 
(HD 150193 is located near the edge of the p Ophiuchus 
cloud). Yudin (2000) estimated the interstellar polarization 
to be 2.5 % at 21 degrees PA in the V- band, and con- 
cluded that the intrinsic optical polarization of HD 150193 
is 4.2%, at 65 degrees PA. By comparing our measured PA 
and Yudin's intrinsic PA, we conclude that most of the po- 
larization we detect in the J band must be intrinsic to the 
disc, with a large amount of polarizing material in the line 
of sight of the star. We note that the integrated polarization 
that we measure around the system of 3.15% is in agree- 
ment with Whittet et al. (1992), given systematic errors of 
0.1 - 0.2%. 

HD 142666: 

This A8Ve star (Dunkin, Barlow & Ryan, 1997b) was cat- 
alogued as a probable Herbig Ae star by Gregorio-Hetem 
et al. (1992), yet no CS disc has been directly observed 
around it. HD 142666 shows evidence of CS activity; a 
double-peaked Ha profile (Dunkin, Barlow & Ryan, 1997b), 
photometric variability, and a large IR excess attributed to 
CS dust (Sylvester et al., 1996). Meeus, Waelkens & Mal- 
fait (1998) interpreted the variable photometric behaviour 
of HD 142666 as evidence of clumps of CS material crossing 
the line of sight of the star. If the CS structure is flattened, 
then the CS disc would be more likely to be orientated edge- 
on. Yudin, Clarke and Smith (1999) measured 0.5% intrinsic 
optical polarization, considerably larger than those of the 
Vega-like 'prototypes' (0.02% and 0.007% for (3 Pictoris and 
Vega respectively; Leroy (1993)). 

Figure 5, left panel, shows our J-band imaging polarime- 
try of HD 142666. HD 142666 presents a polarization pat- 
tern with vectors aligned at ~ 75 degrees PA and a 1.2 % 
mean polarization. As with HD 150193, the PI radial pro- 
file of HD 142666 shows no evidence of extension (Figure 5, 
right panel). Combining catalogued polarization values and 
Hipparcos distances, Yudin et al. (1999) derived a distance- 
polarization law that predicts a 0.3 % contribution at V due 
to interstellar polarization, with a position angle of about 
100 degrees at optical wavelengths. They estimated a V- 
band intrinsic polarization of p ~ 0.5 % PA of 70 degrees. 
The polarization we measure in the J- band is too high to 
be completely interstellar in nature. In addition, we mea- 
sure a polarization PA similar to that of the intrinsic op- 



tical polarization estimated by Yudin et al. (1999). Yudin 
et al. (1999) also noted that the polarization of HD 142666 
probably increases toward IR wavelengths, as is confirmed 
by our J- band measurements. This implies no significant in- 
terstellar contribution to the 1.2% J-band polarization that 
we measure, so we conclude that the polarization detected 
is intrinsic to the system. 

4.2 H- and K- band imaging polarimetry 

H- and K-band imaging polarimetry was obtained for TW 
Hya, HD 150193 and HD 141569 on April 30th 2003. H and 
K photometric values are shown in Table 2, while H- and K- 
band polarimetric measurements are presented in Table 3. 
We note that HD 150193 appears fainter by ~ 0.4 magni- 
tudes in both the H- and K- band with respect to the 2MASS 
measurements. 

TW Hya's H- and K- band imaging polarimetry are shown 
in Figure 6, left and right panels respectively, whilst Fig- 
ure 7 shows the H- and K- band imaging polarimetry of 
HD 150193 (left and right panels, respectively). Alignment 
and subtraction residuals dominate the central parts of the 
images (radius< 0.4 arcseconds). Typical values for instru- 
mental polarization were < 0.4 percent in both H- and K- 
band (as measured from unpolarized standards). As in the 
J- band, HD 141569 did not show any evidence of extended 
polarized structure and is therefore not shown. The near-IR 
spectral dependence of the polarization for our three stars 
targeted at JHK is shown in Figure 8. 

Figure 9, left panel, shows the JHK PI radial profiles of 
TW Hya. Our data show that the radial distribution of the 
JHK PI all follow a similar behaviour between 0.5 and 1.3 
arcseconds from the star, with the polarizing disc being sig- 
nificantly brighter at H than it is at J and K. There is also 
evidence of a decrease in slope in the PI profiles from 1.3 to 
2 arcseconds, as has been previously seen in direct light by 
Krist et al. (2000) and Weinberger at al. (2002). This reduc- 
tion in slope at r > 1.3 arcseconds is clear in both H and K 
but is less evident at J, where a bump is seen at r~ 1.5 arc- 
seconds at very low level (see Figure 1). This is consistent 
with Weinberger et al.'s (2002) total intensity radial profiles, 
in which the change in slope is more evident at H than it is 
at J. 



5 DISCUSSION 

Our polarimetric imaging survey has successfully detected 
the presence of circumstellar polarized emission around four 
of our eleven targets. According to proposed CS evolution- 
ary scenarios, both the IR excess and polarization should 
decrease with time (Malfait, Bogaert & Waelkens, 1998; 
Waelkens, Bogaert & Waters, 1994; Yudin, Clarke & Smith, 
1999); as a consequence of processes such as planet forma- 
tion, the discs should be cleared out and both the near-IR 
excess and the polarizing effects of CS dust should disap- 
pear. Thus, the detection of considerable amounts of po- 
larizing material around TW Hya, HD 150193, HD 142666 
and HD 169142 is consistent with a youthful nature. On 
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the other hand, the non-detection of polarized emission to- 
ward stars presenting low IR-excesses (Lm/Lstar < 0.1) sug- 
gests that these stars are at a more advanced evolution- 
ary state, probably near to the begining of their MS phase. 
In this context we would have expected to detect polarized 
emission around HD 144432 and HD 143006 (which present 
Lm/Lstar values and near-IR excesses even higher than those 
of HD 142666. The most likely reasons for not detecting their 
discs polarimetrically could be that either we did not inte- 
grate deep enough to detect the scattered light against the 
stellar flux or that the disc is symmetrically distributed but 
unresolved, or a combination of both effects. Dominik et al. 
(2003) showed that the near-IR excesses of Class II sources 
(such as HD 144432, HD 142666 and HD 150193) can be well 
modelled if there is a large presence of micron-sized grains 
in the inner 10 AU of their CS discs. It seems reasonable to 
deduce that the near-IR scattering discs around HD 144432, 
HD 143006 and HD 142666 are too small to be resolved by 
UKIRT. In fact, recent VLT mid-IR interferometric observa- 
tions have successfully measured the radius of HD 144432's 
disc to be of 0.014 arcseconds (Leinert et al., 2004), far too 
small to be resolved by our UKIRT observations. 

Apai et al. (2004) proposed that TW Hya's disc Ks-band 
polarized intensity profile and the J- and H-band total in- 
tensity radial profiles have a similar behaviour, concluding 
that the percentage of polarization is nearly constant be- 
tween 0.9 and 1.4 arcseconds from the star. As reported in 
Section 4.2, our observations show that in this region the 
JHK PI radial profiles are indeed very similar. By com- 
bining our PI radial profiles and the total intensity radial 
profiles from HST (Weinberger et al., 2002), we derive TW 
Hya's disc J- and H-band intrinsic percentage polarization 
as a function of radius (Figure 9, right panel). We find Pj 
and Ph arc roughly constant between 0.9 and 1.3 arcsec- 
onds, with Ph being higher than Pj at all radii. In this 
region, the percentage polarization has peaks at 27.6 ±6.1 
at J and 31.3 ± 1.5 at H, both at a radius of 1.2 arcseconds. 
Beyond 1.4 arcseconds Pj falls-off, while Ph rises with in- 
creasing offset from the star. This is a consequence of the 
steeper J-band PI profile, as the Weinberger et al. (2002) to- 
tal intensity profiles for the 1.1 and 1.6 jj,m NICMOS filters 
have the same r~ 2 ' 6 ± 01 fitted slope. However, the NICMOS 
1.6 /im surface brightness deviates from this single power- 
law beyond 1.4 arcseconds, suggesting that our Ph may be 
overestimated beyond this radius. Nevertheless, we can con- 
clude that our observations, in conjunction with direct light 
images from the HST (Weinberger et al., 2002), indicate 
that the disc scattering and polarizing efficiencies peak in 
the H-band, suggesting there is a substantial population of 
sub-micron sized particles present in the disc (Bohren & 
Huffman, 1983). This would argue against larger grains of 
the type suggested by Roberge, Weinberger & Malumuth 
(2005), although this needs to be investigated more thor- 
oughly by future modelling. 

Our imaging polarimetry for HD 169142 confirms the pres- 
ence of an extended CS structure. Both the Q- and U- J- 
band images showed modulation with angle as was already 
noted by Kuhn et al. (2001), suggesting that the disc is 
oriented near to pole-on. The disc-size of 1.2 arcseconds 
we measure (174 AU at HD 169142's distance) is 0.3 arc- 



seconds smaller than the previous polarized-disc extension 
measured by Kuhn et al. (2001). We note that the effec- 
tive exposure time used by Kuhn et al. (2001) was 2.2 times 
larger than the one we used, which may explain the differ- 
ence in HD 169142's measured disc extension between these 
two UKIRT observations. Kuhn et al. (2001) also detected 
a scattering inhomogeneity within the CS cloud, but this 
was not evident in our data. It is interesting to note that 
HD 169142 has been proposed as a Group I source (follow- 
ing the Meeus et al. (2001) classification) and that its SED 
has been successfully modelled with a flaring-disc geometry 
(Dominik et al., 2003) and outer disc-radius of 100 AU. 

Based on the H-band disc-radius of 1.3 arcseconds measured 
coronographically with Subaru by Fukagawa et al. (2003), we 
would have expected to resolve the disc around HD 150193, 
at least in our H- and K-band observations. However we find 
no signs of extension when comparing our PI radial profiles 
of HD 150193 with the total intensity radial profiles of the 
PSF calibrator star. We believe our non-detection of the 
extended disc is due to sensitivity effects, since the Subaru 
observations carried out with a nine times longer integration 
time and a four times larger telescope collecting area. Even 
though we failed to detect the disc reported by Fukagawa 
et al. (2003), we believe there must still be an axisymmet- 
ric structure causing the relatively large linear polarization 
toward this star. The alignment of the polarization vectors 
at 57 degrees PA is likely to be intrinsic, and not associated 
with interstellar polarization (as discussed in Section 4.1). 
To produce the aligned vectors such structure must be (i) op- 
tically thick at J, (ii) inclined to the plane of the sky (i.e. not 
face on) and (iii) too small to resolve with UKIRT. So this 
could be an inner, optically thick component of HD 150193's 
disc. We also note that HD 150193 shows the steepest spec- 
tral dependence of polarization among our three JHK tar- 
geted stars, with its polarization PA being constant within 
the errors at all wavelengths (as can be seen in both Fig- 
ure 7 and in the values presented in Tables 3). Adopting 
the binary PA of 225 degrees measured by Fukagawa et al. 
(2003), we note that there is an apparent alignment between 
polarization PA (57 degrees) and the axis between the two 
stars, although this could be just a coincidence. Neverthe- 
less, this behaviour has already been reported in the V-band 
by Maheswar, Manoj & Bhatt (2002), who suggested that 
it is common to several other HAeBe binary systems. 

We have found for HD 142666 similar evidence of a CS 
scattering structure to that for HD 150193. Along with 
HD 150193, HD 142666 exhibits a large value of linear polar- 
ization compared to the rest of our sample and also presents 
aligned vectors suggesting that the CS disc was unresolved. 
The similarity between both stars in terms of their infrared 
excesses and their polarization makes it plausible to suggest 
a similar structure in their discs. This is consistent with 
the SED modelling carried out by Dominik et al. (2003), in 
which both discs are estimated to be comparable in size and 
mass. 

Despite the limited size of our sample, we note that our 
polarimetric observations detected the presence of extended 
polarizing discs only around systems that have had their 
SEDs fitted with a flared disc geometry (Calvet et al., 2002; 



© 2004 RAS, MNRAS 000, 1-7 



6 Hales et al. 



Dominik et al., 2003). This points toward an observational 
selection effect predicted by Whitney & Hartmann (1992), in 
which Monte Carlo scattered light models show that a flaring 
disc scatters up to two orders of magnitudes more light than 
a flat disc of similar size and mass. Detailed scattered-light 
modelling of our polarimetric data is, however, required to 
constrain the physical parameters of the detected CS discs. 
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Table 1. Programme stars observed 



Name 


Sp class 


Lm/L* 


Distance [pc] 


V 


1 1 1 ^ yoouu 


r\.o ve 


8 A. 1 n — 2 


1 7 


o.ou 


un 1 9^1 fin 


\jO V 


4 4-10 - 3 


1^7 
10. ( 


o.uz 


HD 123356 


G1V 


5-10 -2 


41 


9.7 


HD 141569 


AOVe 


8.4-10~ 3 


99 


7.0 


HD 142666 


A8Ve 


0.34 


114 


8.65 


HD 143006 


G5Ve 


0.56 


82 


10.21 


HD 144432 


A9/F0Vc 


0.48 


200 


8.17 


HD 145263 


FOV 


2-10~ 2 


116 


8.95 


HD 150193 


A1V 


0.37 


150 


8.88 


HD 169142 


A5Ve 


8.8-10- 2 


145 


8.15 


TW Hya 


K7 Ve 


0.25 


56 


11.1 



Spectral classifications and IR excesses are from Sylvester et al. (1996), Dunkin, Barlow & Ryan (1997a) and Sylvester & Mannings 
(2000). An 'e' after the luminosity class indicates the presence of Ha emission. HD 141569, HD 123356, HD 144432, HD 150193 and 
HD 98800 have binary companions at separations of 6.8, 2.2, 1.2, 1.1 and 0.8 arcscconds respectively (Augereau et al., 1999; Sylvester 
& Mannings, 2000; Perez et al., 2004; Mora et al., 2001; Prato et al., 2001). Distances are from the Hipparcos catalogue (Perryman et 

al., 1997). 

Table 2. Photometric measurements 

name Filter Magnitude Magnitude FWHM FWHM 

2MASS reference 
[arcseconds] [arcscconds] 



HD 98800 


J 


6. 


36 


± 





.02 


6 


.39 


± 


0. 


.02 





84 





.87 


HD 123160 


J 


5. 


78 


± 





.03 


5 


.81 


± 


0. 


02 


1 


.30 


1 


.23 


HD 123356 


J 


5. 


94 


± 





02 


5 


.89 


± 


0. 


02 


0. 


,94 





.93 


HD 141569 


J 


6. 


88 


± 





.02 


6 


.87 


± 


0. 


03 





85 


1 


.09 


HD 141569 


H 


7. 


05 


± 


0. 


.02 


6 


.86 


± 


0. 


02 





.41 





.45 


HD 141569 


K 


6. 


84 


± 





02 


6 


.82 


± 


0. 


02 





,38 





,39 


HD 142666 


J 


7. 


32 


± 





.02 


7. 


.35 


± 


0. 


02 





.93 





,92 


HD 143006 


J 


8. 


12 


± 





.03 


8. 


35 


± 


0. 


02 





.86 





.95 


HD 144432 


J 


7. 


13 


± 





.03 


7. 


09 


± 


0. 


02 





85 





.84 


HD 145263 


J 


8. 


10 


± 





02 


8. 


.08 


± 


0. 


02 





.69 





.67 


HD 150193 


J 


6. 


94 


± 





.02 


6 


.94 


± 


0. 


02 





66 





.65 


HD 150193 


H 


6. 


63 


± 





.02 


6 


.21 


± 


0. 


02 





.53 





.51 


HD 150193 


K 


5. 


92 


± 





02 


5 


,47 


± 


0. 


02 





39 





,37 


HD 169142 


J 


7. 


44 


± 





02 


7. 


.31 


± 


0. 


02 


0. 


66 





.58 


TW Hya 


J 


8. 


16 


± 





.02 


8. 


.21 


± 


0. 


02 





82 





.88 


TW Hya 


H 


7. 


.58 


± 





.02 


7. 


.55 


± 


0. 


02 





51 





.46 


TW Hya 


K 


7. 


29 


± 





02 


7. 


29 


± 


0. 


02 





53 





,49 



Photometry was performed using 3.5 arcsccond apertures and calibrated using observed UKIRT standards. Photometric magnitudes 
from the 2MASS point-source catalogue have been included. FWHMs measured both in the targeted stars and repective PSF references 

are presented. 
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Table 3. Polarimctric measurements 



name Filter P [%] 6 



HD 98800 


J 





21 


± 





02 


-70 ± 11 


HD 123160 


J 





12 


± 





01 


80 ± 13 


HD 123356 


J 





26 


± 





01 


84 ± 9 


HD 141569 


J 





55 


± 





02 


-77 ± 9 


HD 141569 


H 





22 


± 





02 


-71 ± 8 


HD 141569 


K 





08 


± 





02 


-59 ± 10 


HD 142666 


J 


1 


32 


± 





02 


75 ± 3 


TTT"\ i ^onru^ 

HD 143006 


J 







± 





02 


6 ± 12 


HD 144432 


J 





49 


± 





02 


3 ± 8 


HD 145263 


J 





37 


± 





03 


23 ± 11 


HD 150193 


J 


3 


14 


± 





02 


57 ± 4 


HD 150193 


H 


2 


04 


± 





01 


52 ± 5 


HD 150193 


K 


1 


41 


± 





03 


48 ± 6 


HD 169142 


J 





22 


± 





03 


-37 ± 6 


TW Hya 


J 





23 


± 





02 


51 ± 6 


TW Hya 


H 





18 


± 





02 


56 ± 7 


TW Hya 


K 





10 


± 





02 


50 ± 8 



Integrated percentage polarizations were computed using the polarimetry tool-box in GAIAt , which allows vectors within a chosen 

aperture to be selected for integration f(http : //www. starlink.rl . ac .uk/). This calculates the mean Stokes Q and U values and 
returns the polarization integrated over the annular region defined by 0.4 <radius< 3.5 arcscconds. This excludes the inner regions 
which arc usually dominated by alignment and subtraction residuals. Polarization position angles are in degrees and measured from 

North to East. 



Radius [AU] 

28 56 84 112 




2.0 1.5 1.0 0.5 0.0 -0.5 -1.0 -1.5 -2.0 ° 05 1 9 

RA Offset (arcsec) 

Radius [ArcsecJ 

Figure 1. Left panel: J-Band imaging polarimetry of TW Hya. Polarization vector maps are superimposed upon total intensity maps, 
showing both contours and grey-scale images. Polarization measurements have been binned in 3 X 3 square bins and a 2a cut applied to 
avoid spurious values. North is up and East is to the left. Right panel: Normalized radial polarized intensity distribution around TW 
Hya (dotted line) compared to the total intensity radial profile of the PSF calibrator star HIP 54690 (solid line). The data have been 
binned over 1-pixcl width rings and the error bars correspond to the statistical errors of each bin. Both x- and y- axis are plotted on 
logarithmic scale. 
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Radius [AU] 
73 



145 




2.0 1.5 1.0 0.5 0.0 -0.5 -1.0 -1.5 -2 
RA Offset (arcsec) 



0.2 



0.5 1 
Radius [Arcsec] 



Figure 2. Left panel: J-band imaging polarimctry of HD 169142, presented in the same format as for Figure 1. The polarization 
vector map of HD 169142 suggests the marginal detection of a centro-symmetric pattern. Right panel: The J-band PI radial profile for 
HD 169142 (dotted line). The solid line corresponds to the total intensity radial profile of the PSF reference star, FS 140. Azimuthal 
values have been binned over 1-pixel width concentric rings. The error bars represent the statistical errors. Both x- and y- axis are plotted 
on logarithmic scale. 





2.0 1.5 1.0 0.5 0.0 -0.5 -1.0 -1.5 -2.0 
RA Offset (arcsec) 




2 4 
Angle [Radian] 



Figure 3. Q- Stokes image of HD 169142 (left panel). HD 169142's Q- Stokes image shows sinusoidal modulation with angle (right 
panel). The crosses mark the data points binned over 18 degree bins, centred on HD 169142, using an annulus extending between 0.6 
and 2.2 arcseconds. The solid curve represents the best sinusoidal fit to the data. 
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Radius [AU] 




RA Offset (arcsec) 

Radius [Arcseconds] 

Figure 4. Left panel: J-Band imaging polarimetry of HD 150193. Polarization vectors are superimposed upon total intensity maps 
shown as both grey-scale and contour images. The polarization measurements have been binned in 3 X 3 square bins and a 2<t cut applied 
to avoid spurious values. North is up and East is to the left. Right panel: HD 150193 A's PI radial profile (dotted line) versus the 
total intensity radial profile of the PSF reference star HIP 80425 (solid line). The region containing HD 150193 B was excluded when 
computing HD 150193 A's PI radial profile. Both x- and y- axis are plotted on logarithmic scale. 



Radius [AU] 




RA Offset (arcsec) 

Radius [Arcseconds] 

Figure 5. Left panel: J-Band imaging polarimetry of HD 142666. North is up and East is to the left. Right panel: HD 142666 's PI 
radial profile (dotted line) versus the total intensity radial profile of the PSF reference star HIP 77815 (solid line). Both x- and y- axis 
are plotted on logarithmic scale. 
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~h [Microns] 



Figure 8. Ncar-IR polarization spectral dependence for HD 150193, HD 141569 and TW Hya (solid, dotted and dashed lines respectively). 
The error bars correspond to the instrumental errors. 
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Radius [AU] 
56 



Mean background level 



J-band 
H-band 
K-band 



Radius [AU] 
56 



J-band 
H-band 



Radius [Arcsec] 



Radius [Arcsec] 



Figure 9. Left panel: J-, H- and K-band polarized intensity radial profiles of TW Hya (dotted, dashed and solid lines respectively). 
Both x- and y- axis are plotted on logarithmic scale. Right panel: J- (circles) and H-band (triangles) intrinsic percentage polarization 
of TW Hya's disc, derived using total intensity coronographic imaging from HST (Weinberger et al. (2002)) and our UKIRT imaging 
polarimetry. The error bars correspond to derived statistical errors. 
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